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Abstract 



General prescriptions for evaluation of coefficients at arbitrary pow- 
ers and logarithms in the asymptotic expansion of Feynman diagrams 
in the Sudakov limit are discussed and illustrated by two-loop exam- 
ples. Peculiarities connected with evaluation of individual terms of the 
expansion, in particular, the introduction of auxiliary analytic regular- 
ization, are characterized. 
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1. The simplest explicit formulae [1-3] (see brief reviews in for the asymptotic 
expansion of Feynman diagrams in various off-shell limits of momenta and masses, 
when the momenta are considered either large or small in the Euclidean sense, have 
been generalized to some typically Minkowskian on-shell limits [5-7], in particular, 
to the Sudakov limit. The prescriptions for these limits have been formulated using 
(pre) subtractions in a certain family of subgraphs of a given graph. 

Recently explicit prescriptions for expanding Feynman integrals near threshold 
have been presented || . This was done with the help of a standard physical strategy 
based on analysis of regions in the space of loop momenta. It should be pointed 
out, however, that this strategy of regions is usually applied only for evaluating and 
summing up the leading logarithms, in particular, in the Sudakov limit (see. e.g. 
||). Note that the information about the leading logarithms is present only in con- 
tributions of some specific regions so that usually one does not consider integration 
in other domains. 

It was argued (and demonstrated for the threshold expansion) in || that it is 
worthwhile to use this strategy for the evaluation of coefficients at any power and 
logarithm for an arbitrary limit. In such extended form, the strategy reduces to the 
following prescriptions: 

(a) consider all the regions of the loop momenta that are typical for the given 
limit and expand, in every region, the integrand in a Taylor series with respect to the 
parameters that are considered small in the given region, 

(b) integrate the integrand expanded, in every region in its own way, over the 
whole integration domain in the loop momenta, 

(c) put to zero any scaleless integral (even if it is not regulated, e.g., within 
dimensional regularization). 

As it was pointed out in ||, it is the step (b) in this procedure that is far from being 
trivial. One may believe that this strategy is legitimate for every limit of momenta and 
masses. For example, it leads to the well-known formulae for asymptotic expansions 
in the case of typically Euclidean limits |], |J (proven in ||) so that we have such 
indirect confirmation at least for them. Note that, for these limits as well as for the 
on-shell limit considered in refs. 0, |6| , the collection of relevant regions is determined 
by subdividing all the loop momenta into large (hard) and small (soft) ones. 

In the present paper, we check, by two-loop examples, this heuristic procedure for 
the evaluation of coefficients at arbitrary powers and logarithms in asymptotic expan- 
sions of Feynman diagrams in the Sudakov limit ]T(J. We shall consider two commonly 
accepted variants of this limit for vertex diagrams with the external momenta p\,P2 
and q = Pi — P2- 

Limit 1. Two external momenta are off shell, p\ = p\ = m 2 = — fi 2 , Q 2 = —q 2 — > 
oo, all internal masses are zero. 

Limit 2. Two external momenta are on shell, p\ = p\ = 0; Q 2 — > oo; some internal 
masses are non-zero. 

We shall calculate the leading power behaviour, including all the logarithms, 



1 



In 7 (q 2 1 m 2 ) , j = 0,1,2,3,4, of the massless planar diagram Fig. la in the first 



limit and compare the obtained result with a known explicit expression |TTJ. Af- 
ter this confirmation we shall then apply the above heuristic prescription to the 
non-planar diagram, Fig. lb, (for which no analytical results are known) in Limit 2 
when mi = . . . = = 0, = itiq = m. We shall also use the second example 
to describe techniques for evaluation of individual terms of the expansion. A natural 
way for evaluation terms with the l/(m 2 ) 2e dependence is introduction of an auxiliary 
analytical regularization. In contrast to the planar diagram in the second limit where 
the poles of the first order in the analytical regularization parameter arise (and can- 
cel in the sum of two contributions) 0, we shall meet, for the non-planar diagram, 
poles up to the second order which are present in five contributions. These poles 
are also cancelled in the sum and we obtain a result which exists within dimensional 
regularization. 

2. The Feynman integral for Fig. la can be written as 

d d kd d l 



FlW ' ro - £) = //(F^ 



I + m 2 ){l 2 - 2p 2 l + m 2 ) 

x \ (i) 

(k 2 - 2 Pl k + m 2 ){k 2 - 2p 2 k + m 2 )k 2 {k -I) 2 ' K J 



We use dimensional regularization with d = 4 — 2e. When presenting our results 
we shall omit iir d / 2 per loop and, when writing down separate contributions through 
expansion in e, we shall also omit exp(— 7^e) per loop (7$ is the Euler constant). 
Let us choose, for convenience, the external momenta as follows: 



2 2 
m _ _, m 



Pi=Pi+ QtPz ,P2=P2+ Qz-Pi , Pi = (Q/2, -Q/2, 0, 0), p 2 = (Q/2, Q/2, 0, 0) (2) 

so that pf = m 2 , p 2 = , 2pip 2 = 2p\p 2 = Q 2 . In the given limit, the following 
regions happen to be typical 0: 



hard (h) 

1- collinear (lc) 

2- collinear (2c) 
ultras oft (us) 



k + ~ Q, k- ~ m 2 j Q , k ~ m . 
k_ ~ Q, k + ~ m 2 j Q , k ~ m . 
k~ m 2 /Q. 



Here k± = k$ ± k\, k=(k 2 , fc 3 ). We mean by k ~ Q, etc. that any component of k^ 
is of order Q. 

One should consider any loop momentum k,l, . . . to be of one of the above types 
and allow for various choices of the loop momenta. (Still it is necessary to avoid 
double counting.) Other types of regions give zero contributions, in particular, when 
one of the loop momenta is soft, i.e. k ~ m. However, if some masses of the diagram 
were non-zero then some soft regions would generate non-zero contributions (that 
would start from a subleading order). 
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In the leading order, 1/Q 4 , we obtain contributions from the following nine regions: 
(h-h), (lc-h), (2c-h), (lc-lc), (2c-2c), (us-h), (us-lc), (us-2c), (us-us). In this list, 
regions for the loop momenta k and / in ([!]) are indicated in the first and the second 
place, respectively. 

The (h-h) region generates terms obtained by Taylor expanding the integrand in 
the expansion parameter, m. In the leading order, this is nothing but the value of 
the massless planar diagram at p 2 = p\ = first evaluated in ref. ||13|| . Although the 
result can be expressed in gamma functions for general e with the help of the method 
of integration by parts (this was first done in |l5j), we present it here, for brevity, 
in expansion in e 



d d kd d l 

{I 2 - 2p 1 l)(l 2 - 2p 2 l){k 2 - 2p 1 k){k 2 - 2p 2 k)k 2 {k - I) 2 



H 5^ 29«3) 3^^_ 

^4e 4 24e 2 6e 32 J (Q 2 ) 2 + 2e 1 ; 

All contributions connected with the ultrasoft regions are easily evaluated in 
gamma functions by use of alpha parameters. In the leading order, we have 



d d kd d l 



/ / >-2p 1 l + m 2 ){-2p 2 l + m 2 )(-2p 1 k + m 2 )(-2p 2 k + m 2 )k 2 (k-l) 2 

r(l-e) 2 r(2e) 2 

- e 2(_ m 2)4 e (Q2)2-2 e > 

d d kd d l 



(5) 



r 0) 

°(us-fc) 



u (us-lc) 



(I 2 - 2p 1 l)(l 2 - 2p 2 l)(-2p 1 k + m 2 )(-2p 2 k + m 2 )k 2 l 2 

' r(i + e )r(i- e )r(e) 2 r(- e ) 2 



T(l -2e)(-m 2 ) 2e (Q 2 ) 2 

d d kd d l 



(6) 



;-2pi/)(/ 2 - 2p 2 l + m 2 )(-2pifc + m 2 )(-2p 2 k + m 2 ) 

1 

x 



■2) 



k 2 (l 2 - (2p 1 l)(2p 2 k)/Q 

r(i- e ) 2 r( e )r(2e)r(- e ) _ (1) 

^(us-2c) ■ K 1 ) 



er(l - 2e)(-m 2 ) 3e (g 2 ) 2 - e ~ ^^ s -2c) 

Using alpha parameters, the rest contributions can be presented, for general e, 
through Mellin-Barnes integrals 

c m = [f 

(lc ~ lc) J J (-2p!/)(/ 2 - 2p 2 l + m 2 )(-2piA;)(A; 2 - 2p 2 k + m 2 )k 2 (k - I) 2 

r(e)r(-e)r(2e) 



r(l + e)(-m 2 ) 2e (g 2 ^ 2 
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1 r+i™ , r(s - 3e)r(s + 1 - 2e)T(s + 1 - e)r(e - s)T(-s) 

X / dS — r — 

T(s + 1 - 3e) 

= 

— °(2c-2c) ' 



2tt J-toc V(s + 1 - 3e) 

a, / / (!''/,•( r'/ 



°(1 



1 L '- / " / ' (I 2 - 2p 1 l){l 2 - 2p 2 l){-2p l k){k 2 - 2p 2 k + m 2 ) 

i r(e)r(-e)r(i - e) 



k 2 (l 2 - (2p 1 k)(2p 2 l)/Q 2 ) T(l - 2e)(-m 2 Y(Q 2 ) 2 +t 

1 /■+*«> r(s + l)r(s - e)r(s + 1 + e)T(-e - s)T(-s) 
x — / ds 




2tt 7-ioo r(s + 1 - 2e) 

= 4% • ( 9 ) 

We imply the standard way of chosing contours: the UV poles are to the right and 
the IR poles to the left of them. The above Mellin-Barnes integrals are expanded in 
e by shifting the contours and picking up residua at points where UV and IR poles 
glue together when e — > 0. As a result we obtain 

°(lc-lc) °(2c-2c) °(lc-/i) °(2c-/i) 

-4C(3)L-^.(10) 

where L = ln(Q 2 /fi 2 ) and we have put \i — 1, for brevity. (Note that, in individual 
contributions, one has both ln(Q 2 //i 2 ) and ln(/^ 2 ).) 

Collecting all nine contributions together we observe that the poles in e which turn 
out to be of very different (UV, IR and collinear) nature cancel, with the following 
result 

(Q 2 ) 2 F l( g,m,0) g -°° \l* + P 2 + 7 -^, (11) 



in agreement with the leading order expansion of the well-known explicit result WL 



3. The expansion of the planar diagram Fig. la, with mi = . . . = = 0, = 
rriQ = m, in Limit 2 was obtained in arbitrary order, following the strategy of sub- 
traction operators, in 0. Note that the same expressions for all contributions of 
the expansion can be obtained with the help of the strategy of regions. The list of 
non-zero contributions, consists, in this language, of (h-h), (lc-h), (2c-h), (lc-lc) and 
(2c-2c) contributions plus a contribution that starts from the next-to-leading order 
and comes from the region where the momentum of the middle line is soft and the 
second loop momentum is considered to be hard. 

Let us now consider the expansion of the non-planar diagram, Fig. lb, in Limit 2. 
The Feynman integral can be written as 



F 2 (Q,m,e) = J J 



d d kd d l 



((k + I) 2 - 2 Pl (k + l))((k + I) 2 - 2p 2 (k + 0) 

X (k 2 - 2 Pl k){l 2 - 2p 2 l){k 2 - m 2 ){l 2 - m 2 ) ' ^ 
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where p\ and P2 satisfy the relations for p\ )2 in the previous section. We shall use 
as well the second choice of the loop momenta when k and I are chosen as momenta 
of lines 3 and 4, respectively, which is obtained by permutation of the masses and 
corresponds to (0) with m\ = m 2 = = rriQ = 0, = WI4 = m. 

Non-zero contributions to the expansion in the leading order are generated by the 
following regions: (h-h), (h-2c), (2c-h), (lc-lc), (2c-2c), (2c-lc), (lc-lc)', (2c-2c)' and 
(us-us)'. As above, we indicate the region for the loop momentum k in the first place 
and for I in the second place. We denote the regions for the second natural choice 
of the loop momenta by prime. The (h-h) contribution is given by the massless 
non-planar diagram. The result, in expansion in e, can be found in fl~3|| : 



f 1 83 C(3) 597r 4 \ 1 

°(h-h) ~ ^4 e2 3e 12Q J (g2)2+2e • W 

The (us-us)' contribution is easily evaluated in gamma functions: 

C (2) = f f 

(us-usy J J (_ 2pi (jfc + l))(-2p 2 (k + l))(-2 Pl k + m 2 ){-2p 2 l + m 2 )k 2 l 2 

= ( g 2)2 i (m2) 4. [r(e)r(2e)r(l - 2e)] 2 . (14) 
The (2c-h) contribution is given by 

d d kd d l 



°(2c-h) 



(I 2 - 2 Pl l + (2p 2 k)(2p 1 l)/Q 2 )(l 2 - 2p 2 (k + l) + (2p 2 k)(2p 1 l)/Q 2 

1 

X (k 2 - 2 Pl k)(l 2 - 2p 2 l)(k 2 - m 2 )l 2 



and the same leading order (h-2c) contribution is obtained by permutation of k and 
I. Using alpha parameters and (twice) Mellin-Barnes representation we obtain 

c m _ c m ( 3 n 2 22 C (3) 16^\ 1 



The (lc-lc) contribution is given by 



(2) 



d d kd d l 



( lc ~ lc ) J J {-2 Pl {k + l))((k + I) 2 - 2p 2 {k + 0) 

X {-2 Pl k){l 2 - 2p 2 l){k 2 - m 2 ){l 2 - m 2 ) ' ^ 

and the (2c-2c) contribution is obtained by permutation of k and I. We should 
also consider similar (lc-lc)' and (2c-2c)' contributions with the second choice of 
the loop momenta. The corresponding expressions are obtained by permutating the 
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masses (see above). The fifth non-zero contribution of the collinear-collinear type 
originates from the (2c- lc) region. It happens that these contributions are regulated 
dimensionally only in the sum. It is convenient to introduce an auxiliary analytic 
regularization into lines 3 and 4 by ^ k 2_ 2pik ^i+xl^2_ 2p2 ^i+x 2 • ^ n con t ras t to the planar 
two-loop diagram in this limit ||, we meet, in this example, poles in Xi up to the 
second order. In particular, the (2c-lc) contribution is evaluated in gamma functions, 
for general e: 



°(2c-lc) 



d d kd d l 



(-2 Pl l + (2p 2 k)(2 Pl l)/Q 2 )(-2p 2 k + (2p 2 k)(2 Pl l)/Q 2 ) 

1 

X (P - 2p 1 k)(l 2 - 2p 2 l){k 2 - m 2 ){l 2 - m?) ' 
T(x 1 )T(x 2 )T(-x 1 - e)T(-x 2 - e)r(x 1 + e)r(x 2 + e) 

r(i + xi)r(i + x 2 )r(- e ) 2 (-m 2 )^+^+ 2e (g 2 ) 2 ' 



;is) 



Using the technique of alpha parameters and Mellin-Barnes representation for other 
four (e-e) contributions, we obtain, for each of them, a result in expansion in X{. Then 
we switch off the analytic regularization (first, x 2 —> X\ and then x% — »• 0), observe 
that, in the sum of all the five contributions, the singular dependence in Xi drops out 
and obtain the following result in expansion in e: 

(n 2^2 \ r {2) , r (2) r (2) r (2) r (2) 

W ) [^(lc-lc) °(2c-2c) °(lc-lc)' °(2c-2c)' °(2c-lc)_ 

_ 19 9 / 2 11tt 2 \ 1 /Stt 2 97C(3)\ 1 
-I?~2? L+ [ L ~^)2T 2 ~{A + ~^)~e 

4 La + *® L - ? W> (19) 

where L = ln(Q 2 /m 2 ) and we have put m — 1, for brevity. 

Collecting all the leading order contributions we see that the poles in e are canceled 
and we arrive at the following result: 

(QrF 2 (Q,m,0) Q -°° T^ 4 -y L2 + 2 °C(3)£-^- (20) 

On the expense of a computer algebra, it is possible to extend this result to any 
order in 1/Q 2 . 

Acknowledgments. V.S. is grateful to M. Beneke, K.G. Chetyrkin and A.I. Davy- 
dychev for useful discussions. 
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Figure 1: (a) Two-loop planar vertex diagram, (b) Two-loop non-planar vertex 
diagram. 
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